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Oh ' Abstract 
< 



In this paper, we will prove the existence of infinitely many positive solutions to 
the following supercritical problem by using the Liapunov-Schmidt reduction 
method and asymptotic analysis: 

f Au + uP + f{x) =0, n > in i?", 
■ 
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1. Introduction and statement of the main results 

The purpose of this paper is to establish the existence of infinitely many positive 
solutions to the following inhomogeneous problem 

r Au + uP + fix) = 0, 



o 
o 



. 5^ , I M > in i?", \im\x\_^^u{x) = 0. 

Where n > 3,p > ^r|,A = J2i=i{d'^ / dx'f) is the Laplacian operator and f{x) G 
C;°'"(i?") with / > everywhere in R^, / ^ 0. 

Inhomogeneous second-order elliptic equations defined in entire space arise natu- 
rally in probability theory in the study of stochastic processes, the eqns. (1-1) in 
particular appeared recently in a paper by Tzong-Yow Lee [9] establishing limit the- 
orems for super-brownian motion. In that paper, existence results for the eqns. (1.1) 
were obtained in the case where the inhomogeneous terms is compactly supported and 
is dominated by a function of the form ^-^^^^^(„.2)p where C{n,p) > is sufficiently 
small. In addition to Lee [9], (1.1) has been studied by Pokhozhaev [11] and Egnell 
and Kay [10]. Pokhozhaev obtained radial solutions when the inhomogeneous term is 
radially symmetric about the origin and satisfies certain integrability conditions. Also, 
Egenll and Kay worked on an equation similar to that of (1.1) but have a small positive 
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parameter £ as a coefficient to inhomogeneous term f{x). In a recent paper [1], exis- 
tence results were obtained for a large class of functions for the inhomogeneous term, 
i.e., 

< fix) < _^ LP 



for all X & R^, where 



2 2 1^— 

4t(''-^-j4i)]'-'. 



Using sub-super solution method, Bae and Ni estabhshed the following infinite 
multiplicity result for the equation 

f Au + uP + fif{x) =0, -u > in R", 
\ lim\^\^^u{x) 0. 

where > is a parameter. 

Theorem A'^l: (i) Let p > Pc- Suppose that near oo 

max(±/(a;),0) < |a;|-''±, 

where g+ > n — A2 and g_ > n — X2 — ^33; ■ Then, there exists /x* > such that for 
every fi G (0,yU*), equation (1.2) possesses infinitely many solutions with asymptotic 
behavior at 00 

(ii) Let p = Pc- Then, the conclusion in (i) holds if we assume in addition that 
either / has a compact support in i?" or / does not change sign in i?", where 



Pc 

and 



{n-2r-4n+4^n^-{n-2f „ ^ in 

(n-2)(n-10) "' ^ 

00 3 < n < 10. 



, (^-2-^) + V(^-2-,4i)^-8(n-2-^4T) 
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Besides, for a more general case (i.e., + \ij is replaced by K{x)uP -|- /x/), similar 
results are obtained, see [2, 8, 12]. However, these results are established require 
p > Pc,n > 11. But the case of n < 10 and < p < Pc is open. Note that in this 
case, the method of sub-super solutions breaks down. In this paper, under reasonable 
conditions on /, we establish that when p > (1.1) has a continuum of solutions. 

The main difficulties in establishing existence of (1.1), in addition to the noncom- 
pactness of the domain and the presence of an inhomogeneous term, are the lack of 
(local) sobolev embedding suitably fit to a weak formulation of this problem. So direct 
tools of the calculus of variation are not appropriate for (1.1). 

Instead of using sub-super sohition method (which limits the applicability on the 
exponent p), we use asymptotic analysis and Liapunov-Schmidt reduction method to 
prove Theorem 1. This is based on the construction of a sufficiently good approximation 
solution. It is well known that the problem 



Aw + wP ^0 in R 



(1.2) 



where p > possesses a positive radially symmetric solution which reduces 

to the equation 

71 — 1 

w" + w' + wP = 0. (1.3) 

r 

This equation can be analyzed through phase plane analysis after a transformation 

2 

introduced by Fowler [9] in 1931: v{s) = rp-^w{r),r = e*, which transforms equation 
(1.3) into the autonomous ODE 

v" + av' - pv + yP ^0 (1.4) 

where 

4 2 , 
a — n — 2 , /3 = (n 



p — 1 p — 1 p— 1 

Since a, j5 are positive for p > the Hamiltonian energy 

E{v) = -v'^ + -^v^^^ - ^v^ 
^ ' 2 p+1 2 

strictly decrease along trajectories. Using this it is easy to see the existence of a 
heteroclinic orbit which connects the equilibria (0,0) and (0,/9~) in the phase plane 
{v, v'), which corresponding respectively to a saddle point and an attractor. A solution 

of (1.4) corresponding this orbit satisfies v{—oo) — 0,i'(+oo) = (3^^ and w{r) — 

2 

r p-iv(logr) solves (1.3) and is bounded at r = 0. Then all radial solutions of (1.3) 
defined in all i?" have the form 

wx{x) := X^w{X\x\), A > 0. 

We denote in what follows by w{x) the unique positive radial solution 

Aw + wP^O in R", w{0) = 1. 

At main order one has 

2 

w{r) Lr p-^ as r — >■ oo, 
which implies that this behavior is actually common to all solutions wx{x). 

2 

By the change of variable A p-^u{j) the (1.1) becomes 

Au + uP + fx{x) = 0, in u> 0, lim u{x) = 0, (1.5) 

where fx{x) — A f{f)- In fact, if / is assumed to satisfy the asymptotic behavior 

f{x) — o{\x\ p-i) as \x\ — >■ +00, 

then we observe that away from the origin /^(a;) — )■ as A — )■ 0. Thus (1.1) may be 
regarded, away from the origin, as small perturbations of problem (1.2) when A > is 
sufficiently small. Prom (1.5), we find that (1.1) "hide" a parameter which indexes a 
continuum of solutions which asymptotically vanish over compact sets. 



Our main result is as follows: 



Theorem 1. Let p > ^^,n > 4. Assume that 

^ n— 3' — 

f{x) = r]{x)fi{x), 

where ri{x) e C^iR""), 0<ri{x) <1, 

ri{x) — for \x\ < Ri,ri{x) — 1 for \x\ > -Ri + 1, 

> fixed large enough) and fi{x) G C°'"(i?"), < \x\~^ with > 2 + 
Then problem (1.1) has a continuum of solutions u\{x) (parameterized by A < Aq, 
where Aq is a fixed number ) such that 

lim u\{x) = 

uniformly in i?" \ {0}. The same results holds when < p < provided that / is 
symmetric with respect to n coordinate axis, namely 

f{xi, ...,Xi, ...,Xn) = f{xi, ... - Xi, ...Xn), for aU i = 1, n. 

The idea is to consider w{x) as an approximation for a solution of (1.1), provided 
that A > is chosen small enough. To this end, wc need to study the solvability of the 
operator A + pw'^'^ in suitable weighted Sobolev space. Recently, this issue has been 
studied in Davila-del Pino-Musso [4] and Davila-Pino-Musso-Wei [5]. 

Throughout the paper, the symbol C denotes always a positive constant indepen- 
dent of A, which could be changed from one line another. Denote A ~ S if and only if 
there exist two positive numbers a, b such that a A < B < bA. 

2. The solvabihty of Hnearized operator A + pw^'^ 

Our main concern in this section is to state the results concerning the existence of 
solution in certain weighted spaces for 

A(/) + pwf- V = /i ini?'*, (2.1) 

where w is the radial solution to (1.2) and his a, known function having a specific decay 
at infinity. We are looking for a solution to (2.1) that is turn out to be a perturbation 
of w, it is rather natural to require that it has a decay at most the same as that of 

2 

w, namely 0{\x\ p-i) as \x\ —> oo. Of course wc would also like be bounded on 
compact sets. As a result, wc shall assume that h behaves like this but with two 

2 o 

powers subtracted, that is, h — 0{\x\ p-^ ) at infinity. 

Now, define some weighted L°° norms as follows (adopted from [4,5]): 

2 

= sup |a;|'^|0(a;)| + sup \x\p-^ \4>{x)\, 

\x\<l \x\>l 

and ^ 

||/i(x)||** = sup |x|^^'^|/i(a;)| + sup |x|^^~|/i(a;)|, 

\x\<l \x\>l 

where a > will be fixed later as needed. 



The following lemmas and remarks on the solvability are due to Davila-del Pina- 
Musso [4] and Davila-del Pino- Musso- Wei [5]: 



Lemma 2.1. Assume that p > ^^,n > 4. For < a < n — 2 there exists a 

^ 71 — 3' — 

constant C > such that for any h with ||/i(a;)||^,* < oo, equation (2.1) has a solution 
4> = T{h) such that T defines a linear map and 

||T(/i)||. <C||/i||... 



For the sake of completeness, we give the main idea of the proof of Lemma 2.1 as 
follows (for the details, see [4, 5]): 

Let ©jfc, /c > be the eigenfunction of the Laplace-Beltrami operator — A^n-i on the 
sphere 5*""^ with eigenvalues repeated according to their multiplicity, normalized 
so that they constitute an orthonormal system in L^(5'""^). We let 6o be a positive 
constant, associated to the eigenvalue and ©j, 1 < i < n is an appropriate multiple 
of which has eigenvalue Xi — n — l,l<i<n. we repeat eigenvalues according to 
their multiplicity and we arrange them in an non-decreasing sequence. We recall that 
the set of eigenvalues is given by {j{n — 2 + j)\j > 0}. We write h as 

h{x) = ^tz+^hke,{9) (2.2) 

and look for a solution to (2.1) in the form 

0(a;) = E^Z+~0fe©fe(e). (2.3) 

Then 

€ + + ipw'-' - ^)<P = h (2.4) 

Equation (2.4) can be solved for each k separately: 

a: If A; = and p > ||/io||** < +oo then (2.4) has a solution 0o which depends 
linearly on hg and satisfies 

||0o||* < C\\ho\U,. (2.5) 
Indeed in this case this solution is defined using the variation of parameters formula 

/r rr 
Z2,ohos"-~^ds - Z2,o{r) ziflh^s'^'^ ds , (2.6) 

where Zi^, Z2fi are two special linearly independent solution to (2.4) with A; = and 
Hq — 0. More precisely, we take Zi^ — rw' -\- -^i^w and 2^2,0 a linearly independent 
solution. Linearization shows that 

n — 2 

^j,ov) = 0{r ~) as r ^ -Foo, j = 1,2, 

while 

-22,0 ~ r^~" near r = 0. 

Using this and (2.6), we can easily get estimate (2.5) 

b: If /c = l,n > 4 and p > < +00, then we have 

||0i||* < C\\hi\U,. (2.7) 



In this case, we have that the positive function zi := —w'{r) > in (0, +00) solves 
(2.4) with k — 1 and hi — 0. Using this, we then define 0i(r) as 

01 (r) = -ziir) r ziisyh^-'^ds T Zi(t)/ii(t)t"-^ (2.8) 
Ji Jo 

Using this formula and by a simple computer, estimate (2.7) is easily obtained. 

c: Let k > 2 and p > If < +00 (2.4) has a unique solution 0^ with 

||0fc||* < +00 and there exists > such that 

||0fc||* < C||/ife||.. (2.9) 

this case is simpler because the operator 

Lk<l> = 0" + + ipw^-' - ^)0 

satisfies the maximum principle in any interval of the form (5, |),5 > 0. Indeed let 
z = —w', so that 2; > in (0, +00) and it is a supersolution, because 

LkZ—— ^ — ^2; < in (0, +00), 

since Xk > 2n for k > n. We construct a supersolution ip of the form 

ip = Ciz + V, v{r) — —. 

Choosing Ci sufficiently large, we can check that 

Lk'^l^ < — cmin(r~'^~^, r~p^~^) in (0, +00). 

Using this, we can easily obtain (2.9). 

The previous construction and (2.5), (2.7) and (2.9) imply that given an integer 
m > 0, if \\h\\ < +00 satisfies hk = \/k > m then there exists a solution to (3.1) 
that depends linearly with respect to h and moreover 

II0II* < C^||/i||** 

where may depend on only m. Then by using a blow up argument, we can show 
that Cfji can be chosen independently of m. For a detailed proof, we refer the interested 
readers to [5,6]. □ 

Remark 2.1. Prom the above proof, we know the operator T in Lemma 2.1 are 
constructed "by hand" by decomposing h and into suns of spherical harmonics where 
the coefficients are radial functions. The nice property is of course that w is radial, the 
problem decouples into an infinite collection of ODEs. 

Remark 2.2. If p > ^^1, linearized operator A + pw^~^ has a kernel, i.e., 
span{^,i = 1,2, ....n}, in general Sobolev space. However, under suitable weighted 
Sobolev space, the linearized operator A + pw^~^ is invertible, i.e., the kernel is 0. 

3. The proof of Theorem 1 



Let p > ^3^, we will prove Theorem 1 in this section. The main idea is to use 
Proposition 2.1 and a contraction mapping principle. 

We look for a solution of (1.5) of the form u — w + cf), which yields the following 
equation for 

Acf> + pwP-'cf>^N{cf>)-Mx), 

where 

N{(f))^-{w + (l)y + wP + pwP-^(f). (3.1) 

Using the operator T defined in Proposition 2.1, we are led to solving the fixed 
point problem 

= T(iV(0)-/,(x)). (3.2) 

Firstly let us estimate || A^(0) ||*^, ^ depending on whether p > 2 or p < 2. 
Case p > 2. In this case, we observe that 

|7V(0)|<C(«;^-V + |0n. 

Let us work with < cr < . Since 

— p-i 

|0(x)| < C|xr^||0||„ for all \x\ > 1, 

and 

w{x) < C(l + \x\)~^, for all x e i?", 
so we have on one hand 

sup \xf+^^wP-''\(f)\'' < C\\(f)f. (3.3) 

On the other hand, 

101 < Clxj-'^ll^ll*, for aU |x| < 1, 

and therefore. We obtain 

sup \x\'+-w{xY-'\<l^\' < \ml sup Ixl'"" = (3.4) 

\x\<l \x\<l 

Prom (3.3) and (3.4) it follows that 

\\w^-'<P%.<C\\<p\\l (3.5) 

To estimate|| 101^11** we compute 

sup|x|2+^|0(x)r<C||0||^. (3.6) 

\x\<l 

Similarly 

sup|a;|2+^|</.r<||0||^ (3.7) 

|x|>l 

Prom (3.6) and (3.7) it follows that 

%*<CU\t (3.8) 



By (3.5) and (3.8) we have 

\\Nm**<cm\i+\\m- m 

Case p <2. In this case |iV(0)| < C|0|p and hence, if < cr < 

sup|x|^+-|0(x)r<C||0||^ (3.10) 

|x|<l 



Similarly 

sup|x|'+F^|0|^'< 11011^ (3.11) 

|x|>l 

From (3.10) and (3.11) it follows that for any 1 < p < 2 and < cr < ^^j, 

||A^(0)||**<C||0||^ (3.12) 

From (3.9) and (3.12) we have 

l|A^(</')ll**<^^(ll</'ll* + ll«- (3.13) 
Now, we estimate \\fx{ 2^)11** as follows: 

sup |x|2+''|/a(x)| = sup \x\'^+''\fx{x)\ < R'^-^-" ^0 as Ri ^ +00, (3.14) 

|a;|<l XRi<\x\<l 

provided that a = 
Similarly, 



sup \xf+^^\fx{x)\ = o{X). (3.15) 

|a:|>l 



2 



So we have, as A — > 0, 



||/A(a;)||,, = sup \x\'+''\h{x)\ + sup \xf^—^\Mx)\ < Rr'^" . (3.16) 

\x\<l \x\>l 

We have already observed that u — w + cj) is a, solution of (1.5) if satisfies equation 
(3.2). Consider the set 

F = {0 e L°°/ II0II, < p} 

where p e (0, 1) is to be chosen (suitably small) and the operator 

h{cl>) = T{N{cP)-f,{x)). 
We now prove that h has a fixed point in F. For e F we have 



mm* < c\\Ni<i>)\u,+c\ 

by (3.9) and (3.10), if cr = Then we have 

\\H<i>)\u,x<c{p'+ff+Rr'-i<P: 

if we choose p is small enough and Ri is large enough. Hence h{F) C F. 



Now we show that ^ is a contraction mapping in F. Let us take 0i, 02 in then 
\\n{(l>i) - h{(l>2)\\. < C\\N{ct>^) - N{(I>2)U. (3.17) 

Write 

7V(0i) - 7V(02) = D^N{4>){4>^ - 02), 

where hes in the segment joining 0i and 02- 
For < 1, 

|xp+'^|7V(0i)-7V(02)| < |xnr»^iV(0)|||0i-02|U 

while, for \x\ > 1, 

|x|'+i^|7V(0,)-Ar(02)| < |xnL'^7V(0)|||0,-02||,. 

Then we have 

||iV(0i) - iV(02)||,, < Csup \x\'\D^N{4>)\\\4>^ - 02||.. (3.18) 

X 

Directly from the definition of A'^, we compute 

D^iV(0) = p[{w + 0)P~i - wP-\ (3.19) 

Thus 

\D^N{4>)\<C{w^-^\4>\ + W-'). 

For all a; we have 

kpw;P-2|0| < C(||0i||, + II02II.) < Cp. (3.20) 

Similarly, for all x 

\x\^w-' < c{m\i-' + ii02iir') < Cff-\ (3.21) 

Estimates (3.13)-(3.15) show that 

sup \x\^\D^N{^)\ < C{p + if-^). (3.22) 

Gathering relations (3.11), (3.12), (3.16) we conclude that ?i is a contraction mapping 
in F, and hence a fixed point in this region indeed exists. So w + 0;^ is solution of 

iS,u + vF + fx{x) = 0, \n u> 0, lim u{x) = 0, (3.23) 

|a;|— >-+oo 

and 

4>x{x) < C for all xeK^X {0}. 

2 

Thus ux{x) — Ap-1 {w{Xx) + (f)x{Xx)) is a continuum of (1.1) and 

lim u\(x) — 

uniformly in \ {0}. This finishes the proof of the theorem 1. □ 

To conclude, in this paper, instead of using sub-super solution method (which limits 
the applicability on the exponent p), we use asymptotic analysis and Liapunov- Schmidt 
reduction method to solve a open problem. 
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